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@ A method for magnetic field gradient eddy current compensation. 



@ A method is provided for characterizing the spurious, 
time-dependent magnetic field gradient response and for 
compensating therefor. In accordance with the method a small 
sample object is placed away from the system isocenter where 
the gradient to be compensated is non-zero. A magnetic field 
gradient pulse is applied to the sample, followed by the appli- 
cation after a variable time of a 90° radio frequency pulse. The 



phase of the resulting free induction decay (FID) signal is mon- 
itored as a function of total time since the end of the gradient 
pulse. These data are related to the gradient impulse response 
and are fit to a set of (e.g.. two or three) exponentials by re- 
gression techniques wherein an initial estimate is improved in 
one embodiment by chi-squared minimization of a Taylor's 
series expansion about the initial estimate point. 
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A METHOD FOR MAGNETIC FIELD 
GRADIENT EDDY CURRENT COMPENSATION 

Backuround of the Invention 

This invention relates to magnetic resonance (MR) techniques. 
More specifically, this invention relates, in the preferred 
embodiment, to a calibration method for MR magnetic field gradient 
eddy current compensation filters. The invention is applicable to 
magnetic resonance imaging, but is not limited thereto. 

The magnetic resonance phenomenon has been utilized in the past 

in high resolution magnetic resonance spectroscopy instruments by 

structural chemists to analyse the structure of chemical 

compositions. More recently, MR has been developed as a medicai 

diagnostic modality having applications in imaging the anatomy, as 

well as in performing in vivo, noninvasive spectroscopic analysis. 

As is now well known, the MR phenomenon can be excited within a 

sample object, such as a human patient, positioned in a homogeneous 

polarizing magnetic field, B , by irradiating the object with 

o 

radio frequency (RF) energy at the Larmor frequency. In medical 
diagnostic applications, this is typically accomplished by 
positioning the patient to be examined in the field of an RF coil 
having a cylindrical geometry, and energizing the RF coil with an RF 
power amplifier. Upon cessation of the RF excitation, the same or 
different RF coil is used to detect the MR signals, frequently in 
the form of spin echoes, emanating from the patient volume lying 
within the field of the RF coil. In the course of a complete MR 
scan, a plurality of MR signals are typically observed. The signals 
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are used to derive MR imaging or spectroscopic information about the 
object studied. 

The application of magnetic resonance to imaging depends upon 
the use of magnetic field gradients to encode spatial information 

5 within the NMR signal. Any departure from the ideal gradient 
behavior can be expected to introduce image distortion. For 
example, imperfect rephasing of the nuclear spins and an attendant 
loss of signal occurs if the gradients are not constant during 
selective time reversal pulse sequences (i.e., use of 180° 

10 time-reversal RF pulses). This effect compounds in later spin 
echoes of multi-echo (Carr-Purcell-Meiboom-Gill) sequences. In 
addition, if the gradient field is not zero when it should be (due, 
e.g., to residual decay from a previous gradient pulse) the 
unintended phase dispersion can result in distorted spectra in 

15 chemical shift imaging (CSI) sequences as well as inaccurate 
spin-spin relaxation time (T ) determination in multi-echo 
sequences. 

Such gradient distortions can arise if the gradient fields 
couple to lossy structures within the magnet, its cryostat (if the 
20 magnet is of superconductive design) or shim coil system, or the RF 
shield used to decouple the gradient coils from the RF coil* The 
spurious response components derive from induction of currents in 
the ambient structures and/or loss of energy to the shim coils and 
are manifested as multi-compartment relaxation behavior. One 
25 observes, for example, an approximately exponential rise and decay 
of gradient fields during and after, respectively, the application 
of a rectangular current pulse to the gradient coil. 

A scheme has been proposed which, in one embodiment, uses an 
analog pre-emphasis filter in the gradient power supply to shape the 
30 current applied to the gradient coil in such a way that the spurious 
gradient components are reduced. The filter may have a number 
(e.g., two) exponential decay components and adjustable 
potentiometers (e.g., four) which must be set during system 
calibration. A measurement technique is used to sample the 
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uncorrected residual gradient field and calculate the potentiometer 
settings from analysis of the data. Other embodiments are 
envisioned wherein the current pulse may be shaped by techniques 
other than pre-earphasis-type filters, such as by generation of the 

5 desired current pulse shape by a computer. 

Accordingly, it is an object of the present invention to provide 
a method for measuring the uncorrected field and to use this 
information to compensate for the gradient distortion. 

It is another object of the invention to provide an improved 

10 method which uses NMR to measure the gradient response and which 
offers improved sensitivity. 

It is a further object of the invention to provide an IRfR method 
to measure the gradient response by monitoring the free-induction 
decay (FID) signal which reflects the entire integrated gradient 

15 history. 

Summary of the Invention 

A method is provided for characterizing the spurious, 

time-dependent magnetic field gradient response and for compensating 

20 therefor. In accordance with the method a small sample object is 

placed away from the system isocenter where the gradient to be 

compensated is non-zero. A magnetic field gradient pulse is applied 

to the sample, followed by the application after a variable time of 
o 

a 90 radio frequency pulse. The phase of the resulting free 
25 induction decay (FID) signal is monitored as a function of total 
time since the end of the gradient pulse. These data are related to 
the gradient impulse response and are fit to a set of (e.g., two or 
three) exponentials by regression techniques wherein an initial 
estimate is improved in one embodiment by c hi -squared minimization 
30 of a Taylor's series expansion about the initial estimate point. 

Brief Description of the Drawings 

The features of the invention believed to be novel are set forth 
with particularity in the appended claims. The invention itself. 
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however, both as to its organization and method of operation, 
together with further objects and advantages thereof, may best be 
understood by reference to the following description taken in 
conjunction with the accompanying drawings in which: 

FIGURE 1 depicts in block diagram form the major components of 
an exemplary MR system; 

FIGURK 2A shows in greater detail the conventional gradient 
amplifier and gradient coil chain; 

FIGURE 2B illustrates a gradient amplifier and coil chain in 
which in accordance with one embodiment of the invention a 
pre-emphasis filter precedes the gradient amplifier; 

FIGURE 3A depicts a sample positioned away from the isocenter of 
an MR system; 

FIGURE 3B depicts an exemplary pulse sequence in accordance with 
the invention for measuring the spurious component of a gradient 
pulse along one axis; 

FIGURE 4A and 4B are similar to FIG. 3B and depict inventive 
pulse sequences suitable for measuring the spurious gradient 
response; 

FIGURE 5 is similar to FIG. 3B t but shows positive and negative 
gradient pulses; 

FIGURE 6 depicts a preferred embodiment of a pulse sequence for 
measuring the position of the sample and for measuring the spurious 
component of a gradient pulse along one axis; 

FIGURE 7 depicts various data curves for performing a regression 
analysis; and 

FIGURE 8 depicts one exemplary embodiment of a pre-emphasis 
filter. 

Detailed Description of the Invention 

Referring first to FIG. 1, there is shown in block diagram form 
the major components of an exemplary MR system. It should be 
understood, however, that this is an exemplary embodiment of the 
apparatus and that in fact the invention is not apparatus 
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dependent. The overall operation of the system is under the control 
of a host computer system generally designated 100 which includes a 
main computer 101 (such as a Data General MV4000). The computer has 
associated therewith an interface 102 through which a plurality of 
computer peripheral devices and other MR system components are 
coupled. Among the computer peripheral devices is a magnetic tape 
drive 104 which may be utilized under the direction of the main 
computer for archiving patient data and images to tape. Processed 
patient data may also be stored in an image disc storage device 
designated 110. An array processor 106 is utilized for 
pre-processing data and image construction. The function of image 
processor 108 is to provide interactive image display manipulation 
such as magnification, image comparison, and gray-scale adjustment. 
The computer system is provided with a means to store raw image data 
(i.e., before image construction) utilizing a disc data storage 
system designated 112. An operator console 116 is also coupled to 
the computer by means of interface 102 and provides the operator 
with the means to input data pertinent to a patient study as well as 
additional data necessary for proper MR system operation, such as 
initiating and terminating scans. The operator console may also be 
used to display images stored on discs or magnetic tape. 

The computer system exercises control over the MR system by 
means of system control 118 and gradient amplifier system 128. The 
computer communicates with system control 118 by means of a digital 
communication network 103 (such as the Ethernet network) in a manner 
well known to those skilled in the art. The system control includes 
several subsystems such as the pulse control module (PCM) 120, a 
radio frequency transceiver 122, a status control module (SMC) 124, 
and the power supplies generally designated 126 necessary to 
energize the components. The PCM utilizes control signals provided 
by computer 101 to generate digital timing and control signals such 
as the current waveforms used for gradient coil excitation, as well 
as RF envelope waveforms utilized in the transceiver for modulating 
RF pulses. The gradient waveforms are applied to the gradient 
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amplifier system 128 generally comprised of G , G and G 

x y z 

amplifiers 130, 132 t and 134, respectively, each utilized to excite 
a corresponding gradient coil in an assembly generally designated 
136 and which is part of a magnet assembly 146. When energized, the 

5 gradient coils generate substantially constant magnetic field 
gradients G , G and G of the magnetic field in the same 
direction as the main magnetic field, wherein the gradients are 
directed in mutually orthogonal I- , Y- and 2-axis directions of a 
Cartesian coordinate system. That is, if the magnetic field 

10 generated by the main magnet (not shorn) is directed in the z 
direction and is termed B , and the total magnetic field in the z 
direction is referred to as B , then G »dB/dx, G * 

Z X Z J 

3B / 3y and G = 3B / 9z, and the magnetic field at any point 
z z z 

(x, y, z) is given by B(x,y,z) » B q + G^X + G y Y + G^Z. The 
15 G^ gradient has no effect on the plane x ■ 0. Similarly, G y and 
,G £ have no effect on the planes y « 0 and z * 0, respectively. 
The point (0,0,0) is termed "isocenter" and is the point in space 
where none of the gradients have any effect. Isocenter is normally 
situated substantially at the center of the static magnetic field 
20 volume. 

The gradient magnetic fields are utilized in combination with 
radio frequency pulses generated by transceiver 122, RF amp 128 and 
RF coil 138 to encode spatial information into the MR signals 
emanating from the region of the patient being studied. Waveforms 

25 and control signals provided by the pulse control module are 
utilized by transceiver subsystem 122 for RF carrier modulation and 
control of the operating mode; that is, the transmit or receive 
mode. In the transmit mode, the transmitter provides a radio 
frequency carrier waveform modulated in accordance with the control 

30 signals to an RF power amplifier 123 which then energizes RF coils 
138 which are situated within main magnet assembly 146. The NMR 
signals radiated by the excited nuclei are sensed by the same or a 
different RF coil than is used for transmitting. The signals are 
detected, amplified, demodulated, filtered, and digitized in the 
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receiver section of the transceiver. The processed eignals are 
transmitted to the main computer for processing by means of a 
dedicated, unidirectional, high-speed digital link 105 which links 
interface 102 and transceiver 122. 

5 The PCM and SCM are independent subsystems both of which 

communicate with main computer 101, peripheral systems, such as 
patient positioning system 152, as well as to one another by means 
of link 103. The PCM and SCM are each comprised of a 16-bit 
microprocessor (such as an Intel 8086) for processing commands from 

10 the main computer. The SCM includes means for acquiring information 
regarding patient cradle (not shown) position and position of the 
moveable patient alignment light fan beam (not shown). This 
information is used by the main computer to modify image display and 
reconstruction parameters. The SCM also initiates functions such as 

15 actuation of the patient transport and alignment systems. 

, The gradient coil assembly 136 and the RF transmit and receiver 
coils 138 are mounted within the bore of the magnet utilized to 
produce the polarizing magnetic field. The magnet forms a part of 
the main magnet assembly which includes the patient alignment system 

20 148, a shim coil power supply 140, and a main magnet power supply 
142. The shim power supply is utilized to energize shim coils 
associated with the main magnet and which are used to correct 
inhomogeneities in the polarizing magnetic field. In the case of a 
resistive magnet, main magnet power supply 142 is utilized to 

25 continuously energize the magnet. In the case of a superconductive 
magnet, the power supply is utilized to bring the magnet to the 
proper operating field and then is disconnected. In the case of a 
permanent magnet, power supply 142 would not be needed. The patient 
alignment system 148 operates in combination with a patient cradle 

30 and transport system 150 and patient positioning system 152. To 
minimize interference from external sources, the MR system 
components comprised of the main magnet assembly, the gradient coil 
assembly, and the RF transmit and receiver coils, as well as the 
patient-handling devices, are enclosed in an RF-shielded room 
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generally designated 144. The shielding is generally provided by a 
copper or aluminum screen network which encloses the entire room. 
The screen network serves to contain the RF signals generated by the 
system, while shielding the system from RF signals generated outside 

5 the room. A bidirectional attenuation of approximately 100 db. is 
typical in the frequency range of operation. 

Referring now to FIG. 2A, there is shown gradient amplifier 130 
(also shown in FIG. 1) for energizing a gradient coil 12 to produce 
the G x gradient. Under ideal operating conditions, a rectangular 

10 current pulse 14 applied to amplifier 130 would be amplified, as 
indicated at 16, and when used to energize the gradient coil would 
result in a substantially rectangular gradient magnetic field pulse 
18 being produced. However, as described hereinbefore, because of 
the coupling to lossy structures and spurious components due to 

15 induction of eddy currents, in practice, the resulting magnetic 
field gradient has a finite rise time as indicated at 20, and a 
finite decay as indicated at 22. As alluded to hereinbefore, such 
gradient distortions can lead to loss of signal and unintended phase 
dispersion. 

20 In accordance with one embodiment of the invention the gradient 

distortions can be reduced by application of a current pulse 14, as 
shown in FIG. 2B, to a pre-exnphasis filter 24 to predistort the 
current pulse, as indicated at 26, such that the amplified current 
pulse 28 applied to the gradient coil would produce the desired 

25 rectangular gradient pulse 18. Since in typical MR applications, 
gradient pulses are applied in at least each of the axes of the 
Cartesian coordinate system, an MR system for practicing the 
invention would have means functionally similar to that depicted in 
FIG. 2B to achieve correction along all axes. However, for 

30 simplicity the following discussion will refer to correction along 
only one of the axes. 

In order to determine how current pulse 14 in FIG. 2B should be 
predistorted and therefore how pre-emphasis filter 24 should be 
calibrated to achieve the desired distortion, the nature of the 
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gradient distortion to be eliminated must first be measured .and 
analyzed. The manner in which this is accomplished using an KB 
system will be described next. 

FIGURE 3A depicts a sample 30 composed of an HR-active 

5 substance, conveniently a small bottle containing approximately one 
cc of water, positioned at a distance ,- x M from the gradient origin 
designated M 0" (that is, the system isocenter) where the gradient is 
non-rero. Sample 30 will serve as the source of the VMR signals 
used to study gradient distortion using the pulse sequence 

10 diagrammed in PIC. 3B. F I CURB 3B shows only a single gradient 
applied, for example, along the x-axis of an object. As previously 
indicated, the process to be described would in practice be repeated 
for other axes along which gradients requiring compensation are 
applied. 

15 It may be helpful to initially consider an intuitive explanation 

of the principle underlying the present measurement method. As is 

# 

shown in FIG. 3B, a gradient pulse is followed by the application of 
o 

a 90 RF pulse which generates a free induction decay (FID) 
signal. If there are no spurious gradient components, then the 

20 sample object will be immersed in a constant magnetic field for all 
time during the FID. As a result, the instantaneous frequency of 
the FID will be constant as a function of time. If, however, the 
gradient pulse is followed by spurious components as shown by the 
dished curve in FIG. 3B, then the magnetic field at the sample 

25 object will vary during the FID and so will the instantaneous 
frequency of the FID signal. Referring again to FIG. 3B, there 

is shown a gradient pulse being applied just before a 90° RF 
pulse. The resulting free induction decay signal will contain phase 
information whose instantaneous value at time **t M depends on the 

30 time integral of the gradient since the application of the RF 
pulse. 

The FID signal, S(t) which is obtained from a sample with one 
dimensional density distribution P(x) (i.e., integrated over the y 
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and z directions) and spin-spin relaxation time T 0 (x) is: 




5 

where y is the gyromagnetic ratio, 6 (t) represents the gradient 
response of interest and u is a constant resonance frequency 
offset which is discussed in more detail below. If the object is 
small (the effect of object size is discussed fully below), the 
10 phase of the signal S(t) is 

t 

:(t) = yx J G (t')dt' + ie^t • (2) 

o * 

15 It is seen that the phase of the FID signal contains information 
about the gradient response C^Ct). The present invention utilizes 
this information to characterize the spurious gradient response 
using this measured phase. A further goal is to use the measured 
information to pre-corapensate the gradient amplifier's current 

20 waveforms so as to produce the desired gradient fields. Thus in 
principle one need merely extract the phase from the FID signal as a 
function of time in order to characterize the spurious gradient. In 
practice, however, phase dispersion due to both constant B 

o 

inhomogeneity as well as from the gradient-induced inhomogeneity can 
25 lead to a T 2 * decay much shorter than the characteristic 
relaxation times of the gradient. In this case the FID signal will 
decay before the full spurious component of the C x <t) curve has 
been sampled. 

As a result of the Tg* decay, a single experiment such as is 

30 shown in FIG. 3B can only characterize the gradient field for a time 

o 

comparable to T^*. Suppose the 90 RF pulse follows the 
gradient pulse after essentially no delay, and the FID is measured 
for a time T^ msec. (T^ <. T 2 *) immediately following the RF 
pulse. This measurement will characterize the gradient field for 
35 the first ? msec, after the current pulse is shut off. 
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How suppose the application of the RF pulse is delayed by 
msec, after the gradient current pulse is shut off and again sample 
the FID for msec, after the RF pulse. If < T it 

will be possible to combine the data from this and the previous 

5 (T^ » 0) experiment to allow the gradient field to be 
characterized for a longer period of time. This process can be 
repeated so that the full spurious response can be measured. 

The sequence of FIG. 3B is therefore modified as shown in FIG. 
4A so that a variable time delay T 3 is inserted between the 

10 gradient and RF pulse. A number of short acquisitions of length 
<L T 2 * are then obtained with T 3 incremented by A 
T^ between acquisitions. In this way, the entire spurious 
gradient decay can be swept through the A/D acquisition window in 
small overlapping segments. In one implementation, T^ was 

15 selected to be 10 msec, while the change in the delay between 
segments was 8 msec, as shown in FIG. 4B, in which A/D^ t A/D^, 
etc., indicate adjacent, overlapping segments in sampling the 
gradient spurious response curve. The total phase curve 22 is then 
obtained by combining the separate acquisitions 34, 35 and 36 shown 

20 in FIG. 4B. 

In the analysis of the phase data from each acquisition and in 
the process of combining the data from the various acquisitions, use 
is made of the fact that $ (t) (Equation (2)) is continuous. Thus, 
even though the phase at any single point in time can only be 

25 calculated modulo 2 w t values of <><t) much greater than 2* can be 
measured by "unwrapping** the measured phase. Similarly, the various 
segments are combined by forcing 4>(t) to be continuous at the 
segment boundaries. In order to better match the segments at the 
boundaries, it is advantageous to overlap the segments somewhat. 

30 In this way, $(t) can be measured over time spans much greater 

than T^* which are long enough to fully characterire the spurious 
gradient decay. A suitable analysis of these data then provides the 
desired characterization in the form. for example, of a 
multi-exponential decomposition. The ensuing discussion will 
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describe the theory and regression analysis used to obtain the 
decomposition more fully. 

Let the impulse response of the gradient be h(t). Because of 
the ease with which measured parameters can be converted to circuit 
5 component values, and because it is found to be a good 
approximation, it is convenient to approximate the impulse response 
of the gradient system by the sum of a finite number of exponential 
terms. Thus, h(t) is defined by: 

10 h(tj = i(t) + Z H gV*** . (3) 

k=l k 

in which 5 (t) is the Dirac delta function and g^ and X are 
the amplitude and rate constant, respectively, of the k r th 
15 exponential component. The gradient response is obtained by 
convolving the impulse response with the excitation function. For a 
rectangular excitation and pulse the definitions for the time 
parameters shown in FIG. 3A one obtains: 

G(t=0) • G 0 £ a k (l-e" XkT °)-= % G fc _ (<) 



25 where 



**-h'\ , (5 , 



is a dimensionless amplitude coefficient* Equation (4) expresses 
the gradient amplitude just after the end of the gradient pulse with 
30 the assumption that an infinite time has elapsed since any prior 
gradient excitation. After the excitation is terminated. 

G x (t ' ) = G o r k G k e " V • (6) 
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Substituting Equation (6) into Equation (2) yields: 

5 

where T fc » A* 1 is the relaxation time. Thus, if the 
position of the sample x (FIG. 3A) and the offset Au are known, 
and 9 (t) is obtained from measurements of S(t), the coefficients 
ot fc and needed to calibrate pre-emphasis filter 24 (FIG. 2B) 

10 can be obtained by chi-squared minimization techniques as discussed 
below. 

The offset phase A u t derives from frequency synthesizer offset 
and/or field inhomogeneity . In the present embodiment. the 
frequency £ a is determined by making a second measurement with the 

15 sign of the gradient G (t) reversed relative to the preceding 
G x (t) gradient pulse, as shown in FIG. 5. The resulting signals 
S*(t) and S~(t) have opposite phases as suggested by the + and - 
superscripts. Since the phase due to the offset in Equation (7) 
does not change sign when the gradient sign is reversed, the offset 

20 frequency component is removed when the two measurements are 
subtracted. This is analogous to the method disclosed and claimed 
in commonly assigned U.S. Patent No. 4,443,760. 

The manner in which the position "x M of sample 30 (FIG. 3A) is 
measured using KB will be described next. The position of the 

25 sample is determined from the spectrum of a gradient-ref ocused 
spin-echo, as depicted in FIG. 6. 

The pulse sequence of FIG. 6 is similar to that of FIG. 5 except 
that it includes an imaging experiment used to measure M x". The 
experiment includes the initial application of a 90° RF excitation 

30 pulse, followed by application of opposite polarity gradient pulses 
40 and 42 which produce the spin-echo signal. Gradient pulse 42 is 
the readout gradient pulse (wherein » 120 mg/cm) , for example. 
If, as discussed hereinbefore, gradient correction is desired along 
other axes (e.g., y and z), a similar imaging experiment along those 
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axes would be needed to measure the position ("y M and M z") of the 
sample. 

Referring again to FIG. 6. let C p be the amplitude of the 
readout gradient. Then the frequency for the sample at position M x" 
5 is 

o P (8) 

where fa is a constant offset. The offset frequency is removed by 
10 gathering half the views with G p negated and subtracting the 
frequencies of the two half-scans, as described hereinabove. The 
frequencies are found by a standard thresholded moment analysis of 
the spectra. • 

The reason a spin echo is used to measure "x" instead of 

15 geometric means is that (1) the origin (or null point) of the 

gradient may be difficult to find, and (2) the gradient may not be 

perfectly linear, especially at the edges of the field of view. 

Thus, by using magnetic resonance signals to calculate T, any 

errors in G (FIG. 6) will be compensated by the same errors in 
o 

20 G . Only the ratio G /G is important. 
P P o 

Hext, the effect of using an object of finite size is 

discussed. If the sample is not infinitesimal, as assumed in the 
discussion heretofore, the spurious gradient response will cause 
dephasing across the width of sample 30. This effect is like T 2 * 
25 in that it limits the FID duration and can alter the relationship 
between <p (t) and the phase of S(t). 

For simplicity assume a rectangular sample of width 2a centered 
at V. Then the FID is 



30 



t 

x+a iyx S G(t')dt' 
S(t) =1 e 0 dx 

2a x-a . < 9) 
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or 



S(t) - e i0(t) sinc[*(t)a] 

x 



(10) 



Thus, Sit) is amplitude-modulated by the sine function. Mote from 
Equation (10) that the finiteness of the sample does not distort the 
phase measurement 6ince ars (S(t) » $ (t). It can be shown that 
this relationship holds for any symmetric object. Thus, no error is 
made in the measurement if the sample is symmetric even if it is not 
infinitesimal. 

The modulation effect may be estimated by letting the gradient 
during the FID be constant. Then the first null in the amplitude 
occurs at time T^*, where 

XfG o a V = 7 • 



Thus an upper limit to the gradient C is obtained by setting the 

o 

FID duration T^* equal to the A/D window duration. For example, 
if a ■ 8% and if a T 2 * 2 10 msec, is desired, then G o <. 200 
mgauss/cm. 

The final data determination after phase offset correction and 
"x" determination is 



M(t) - ilH i . La k (l.e" T ° Xk ) Tlf (l.e" a 0 

«o G o " (12> 

Thus, coefficients A^, t fc are sought In the expansion 

M(t i } " VSc 0" e tl/T ^ • 1 = 1 N m • (13> 
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where there are H measurements at times t. , and 
m l 



k = a k V- e J T i 



(14) 



5 In one imaging system it is found that two or three exponentials fit 
the data quite well. In other applications, fewer or more 
exponential components may be needed. The techniques used to decide 
how many terms are needed are well known in the art as disclosed, 
for example, in a text by Bevington identified hereinafter. 

10 To characterize the spurious gradient response, the coefficients 

(A^) and ( x fc ) must be calculated. A reasonable approach is to 
perform a least square fit (or minimization). As can be seen 

in Equation' (13), the expression to be fit is not linear in these 
coefficients. Methods for non-linear minimization are well 

15 known. Several are described in the Bevington text. A description 
of the method used to fit data from one particular imaging system 
will now be described. Either two or three exponential terms may be 
fit. 

First, a preliminary fit for the longest time constant term is 
20 calculated. For this preliminary fit, data for late times is used, 
as shown in FIG. 7A, so that the shorter time constant terms can be 
ignored. For time "f greater than the shorter time constants 

-t/T M 

M(t)^A Q -A N e . (15) 



where H is th* exponential term with the longest time constant and 



N 

A 0 = ^ A k . (16 , 



2 

30 A conventional iterative X minimization method is used to 
calculate A q , ^ and x g (Equation (15). The coefficients 
Ag and T g comprise the preliminary fit to the longest time 
constant term. Hext, this longest term is subtracted from the 
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measured data 

"•(t) = M, t ).A N (,.e- t/T »)= I "- , A k 6.e- t/T 0 • <n> 

5 

and this modified data is used to perform a preliminary fit for the 
term with the second longest time constant. Again, a window of data 
is chosen so that the effects of any remaining term can be ignored. 
This is graphically shown in FIG. 7B. Using a procedure similar to 

10 the one described above, A^ ^ and T ^ are calculated. If 
three time constants are being considered, the process of subtractng 
the preliminary fit and computing the coefficients for the next 
shorter term are repeated, as shown in FIG. 7C. 

In this way, preliminary values for the coefficients are 

15 calculated. It is found that further processing, specifically to 

simultaneously fit to all the exponential terms, can improve the 

fit. One simple procedure is to keep the time constants ( T ) 

fixed to their previous estimated values, but to calculate the 

2 

values of the amplitudes (A^) which minimize X ■ As can De 
20 seen in Equation (13), this fit is linear in the amplitude 
coefficients and so it involves a straight forward linear 
least-square fit and the solution of an tfxtf system of equations. 

Another calculational method that has been found to be useful is 
to improve the fit by a Marquardt procedure with a Taylor's 
25 expansion of the x function about the initial guess point 
retaining only terms up to first derivatives, as described in P. R. 
Bevington, "Data Reduction and Error Analysis for the Physical 
Sciences," McCraw Hill Book Co., II. Y. , 1969, pp. 232-235. This 
process works well if the initial guess is reasonably close. A 
30 2Vx2M linear system of equations is solved by standard numerical 
matrix inversion techniques. It is found that digital computational 
accuracy, even when double precision floating point arithmetic is 
used, can limit the accuracy of the fit coefficients. To improve 
the accuracy of the fit, only the time constants resulting from this 
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10 



15 



20 



25 



Marquardt procedure are retained. The amplitude coefficients are 

calculated using the linear fit procedure described above. 

Thus, the complete fit procedure consists of: (a) preliminary 

fit to each of the exponential terms by subtracting the previous 

preliminary fits and using only a subset of the data; (b) retain 

only the time constant coefficients ( T fc ) and calculate new 

amplitude coefficients (A.) using a linear fitting procedure; (c) 
2 

expand the x expression using a Taylor series about the previous 

2 

fit results and calculate a new fit to minimize X ; and (d) 
retain only the time constant coefficients ( T fc ) and calculate new 
amplitude coefficients (A^) using a linear fitting procedure. 

The steps (c) and (d) are repeated until X 2 no longer 
decreases significantly. 

Equation (7) described above, assumes that the gradient response 
has fully recovered since any previous pulse by the time the 
sequence is played out. In practice, this is not the case for 
components with long-time constants T ^ (i.e., comparable or 

larger than the separation time between gradients), and a correction 
must be made as will be described next. 

Suppose fc(t) is modelled with exponential components (H»3) and 
that only the third component (k«3) has a long enough time constant 
for the present effect to be operative. Because of the application 
of the positive, G x (t)f ud negative, G~(t), gradients 
described above, there are two gradients of alternating signs in the 
pulse sequence, as depicted in PIG. 5. For this sequence, the field 
following the positive gradient after dynamic equilibrium has been 
set up is 



where G (t) is the gradient in the absence of overlap of spurious 
o 

components (Equation 7), and, again, only the third exponential 
component of the impulse response is presumed long enough for this 




(18) 
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overlap effect to be important. A similar equation holds for the 
spurious gradient following the negative pulse c"(t). Since 
the total offset-corrected phase is proportional to the difference 
between and c ^( fc >» the total phase due to just the 

5 long component is (from Equation (18) 

where t> 3Q is the infinite TR phase given by Equation (7) for 
k-3. Since, as will be shown, o- 3 is proportional to <t> ^, the 
corrected a 3 becomes 

15 a 3 = a 3 (meas) ( 1 + e"* 3 ^- j e"* 3 ^V* 3 * 6 ] . (20) 

1 - e" 3 R 



where a 3 (meas) is the uncorrected measurement. 

20 

Once the coefficients of the impulse response are determined* as 
discussed above, the deconvolution (pre -emphasis) filter may be 
determined as described next. 

The significance of the impulse response given in Equation (3) 
25 is that if a deconvolution filter *(t) - h(t)" 1 is first applied 
to the excitation, the net response will follow the desired 
excitation. Thus the pre-emphasis filter is obtained as the inverse 
of the impulse response. To calculate this, let 

30 h(t) * 6(t) «(t) (21) 

define the error function c(t) . Then, it may be verified (h*k*6(t)) 
that the appropriate filter is 

35 k(t) « £(t) - c(t) ♦ c* c(t) - e*£*e+... (22) 
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where denotes convolution. Hote that if e is of order a % 5%, 
then the higher order terms e* e, etc., are negligible. It has 
been chosen that only the first term is retained since then the 
pre-emphasis filter, in one embodiment, is just a series of first 
order RC networks with amplitudes and time constants x^. 

One exemplary embodiment of a p re-emphasis filter, such as 
pre-emphasis filter 24, FIG. 2, is illustrated in circuit schematic 
form in FIG. 8. The amplifier is made up of an operational 
amplifier 80 and filter components 82, 83, and 84. The predistorted 
output which is applied to gradient amplifier (e.g., 130, FIG. 2B) 
for energizing a gradient coil. At its input the operational 
amplifier is energized by a current pulse 14 through an input 
resistor Ro and one or more parallel-connected RC networks 82-84. 
The number of RC networks depends on the desired precision in 
pre-emphasis to be applied to the current input 14 to achieve 
compensation in the gradient pulse. Dash lines 86a and 86b suggest 
that additional RC networks may be added as needed. Each RC network 
is comprised at its input of a capacitor (e.g., CI) connected in 
series with an output variable resistor (e.g., Rl). The common 
point between the capacitor and variable resistor is coupled to 
ground through a second variable resistor (e.g., R2). The circuit 
shown in FIG. 8 provides three exponential correction (since there 
are three RC networks). In some applications it has been found that 
two RC networks provide adequate correction (pre-emphasis). 

After the coefficients are calculated, the resistance settings 
for the variable resistors (e.g., R1-R6) in the RC networks can be 
obtained. In principle, if the response of amplifier 24 has been 
perfectly modeled by the circuit calculation, the measured response 
will be perfectly compensated. In practice, of course, the 
component tolerances of capacitors may result in some departure from 
the ideal. In this case an iterative correction step may be used 
wherein the measured response after correction is fit to a new set 
of exponentials using the previous time constants, and the resulting 
error amplitudes for each component are added to the previous 
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values. This iterative process may be repeated, although it has 
been found that in some cases no iterations are needed. 

Although one embodiment of the invention ha* been described 
above which uses an exponential pre-emphasis filter, it will be 

5 recognized that other filter implementations (such as active filter 
designs) may be devised by those skilled in the art. Additionally, 
embodiments are envisioned in which following measurement of the 
decaying spurious component pre-emphasis (or predistortion) of the 
signal energizing the gradient coil may be achieved using software 

10 techniques without the aid of pre-emphasis filters. 

In accordance with such software techniques, for example, once 
k(t) is found using Equation (22), it may be applied by numerical 
convolution techniques to the data used to generate the gradient 
waveform, prior to applying the gradient waveform to the gradient 

15 amplifiers. 

Thus, let w(t fc ) be the desired gradient waveform for a 
gradient to be applied along particular axes. Then the corrected 
waveform which must be applied to the gradient amplifier for the 
actual response to be corrected is 

20 N 

W c (t k ) = Z w(t ) . k(t ) f (23) 

where H denotes the number of points in the sampled filter kernel 

k. This operation may either be performed by a computer (such as 
25 computer 101, FIG. 1) before the waveform is stored in a readout 

memory, or it may be performed in nearly real time by computation 

from a stored wave w(t^). 

While this invention has been described with reference to 

particular embodiments anc examples, other modifications and 
30 variations will occur to those skilled in the art in view of the 

above teachings. Accordingly, it should be understood that within 

the scope of the appended claims the invention may be practiced 

otherwise than is specifically described. 
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CLAIMS 

The invention claimed is: 



1. A method for magnetic field gradient eddy current 
compensation comprising the steps of: 

(a) positioning an BJMR-active sample away from the 
isocenter of an MR scanner apparatus; 

(b) applying along at least one dimensional axis of said 
sample a gradient pulse; 

(c) irradiating said sample with an RP excitation pulse 
within a predetermined time following application of said gradient 
pulse so as to generate an HMR signal; 

(d) measuring the phase of said NMR signal as a function of 
the total time since the application of said gradient; 

(e) relating the measured phase to the impulse response of 
said gradient; and 

(f) compensating for magnetic field gradient eddy current 
using said impulse response. 
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